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ABSTRACT 


A model is presented for estimating the dimensions of delaminations in continuous 
fiber reinforced composite plates subjected to non-penetrating (“low velocity”) impact of 
a solid object. The model is based on dimensional analysis, and results in two simple, 
closed form expressions, one providing the delamination length, the other the delamination 
width. Comparisons of delamination lengths and widths calculated by these expressions 
with experimental data suggest that the model predicts these dimensions with reasonable 
accuracy. 
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I. INTRODUCTION 


Delaminations in composites caused by the impact of a solid object are of great practical 
concern. Delaminations may reduce significantly the strength and stiffness of the material, 
especially when the impacted laminate is under compression. Therefore, in order to utilize 
composites to their full advantage, we must be able to assess their susceptibility to impact. 
Specifically, it is important to know the locations and dimensions of delaminations which 
may occur due to impact. 

Recently, Wu and Springer [1] described a method for estimating the locations and sizes 
of delaminations inside fiber reinforced organic matrix plates impacted by a solid object. Wu 
and Springer’s analysis utilizes a three dimensional, transient, finite element stress analysis 
coupled with a failure criterion based on dimensional analysis. In addition to the delami- 
nation sizes and locations, this analysis also provides detailed behavior of the impactor and 
the impacted plate (e.g., displacement and velocity of the impactor, contact force as a func- 
tion of time, and displacements, velocities, stresses and strains in the plate as functions of 
position and time). The disadvantage of the Wu-Springer analysis is that it requires long 
computation times on mainframe computers. 

Frequently, in practice, the full history of the impactor and the plate is not needed; 
an estimate of the locations and sizes of the delaminations suffices. This investigation was 
undertaken to develop a method which provides this information simply and easily, without 
recourse to very long computer calculations. The present method is based on the failure 
criteria proposed by Wu and Springer, but does not require the time consuming finite element 
calculations. Instead, the method described in this paper results in closed form expressions, 
which readily yield the locations and sizes of delaminations caused by non-penetrating (“low 
velocity”) impact. 
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II. THE PROBLEM 


We consider a rectangular plate (thickness h) made of a continuous fiber-reinforced, 
organic matrix composite (Figure 1). The plate consists of n plies or lamina. The ply 
orientations are arbitrary and need not be symmetric with respect to the mid-surface of the 
plate. Perfect bonding between each ply is assumed. 

The plate is supported along all four edges. Each edge may be either clamped or simply 
supported. A solid object of mass m, speed V, and “nose” radius r impacts the center of 
the plate (Figure 1). It is assumed that during impact the impactor does not penetrate the 
surface. The objective is to develop a method which can be used to estimate the locations 
and sizes of delaminations caused by the non-penetrating impact. 


III. DELAMINATION 

The shape of the delamination is irregular. However, previous data have shown the 
major axis of the delamination is, in general, parallel to the fiber direction of the lower ply 
[2]. Therefore, we characterize the delamination between any two plies by the length tj) and 
width wd which are taken to be the dimensions of the delamination in the directions parallel 
and perpendicular to the fiber direction in the lower ply, i.e., in the ply further from the 
impactor (Figure 2). 

By using dimensional analysis, Wu and Springer [1] obtained the following expressions 
for the delamination lengths and widths 



Kj c and to axe ply properties which must be determined by tests. Kj c is the Mode 
I critical stress intensity factor (opening mode), to is the critical crack length given by 
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transverse tensile strength measurements [1, 3] 


2Kl 

irSj^rp 


where Stt is the transverse ply tensile strength. 


( 3 ) 


Q* and Dgg are laminate properties which can be calculated from laminated plate theory. 
Q* is proportional to the difference between the stiffnesses of the two plies adjacent to the 
delamination in the direction of the fibers of the “lower” ply (Figure 2). The parameter Q * 
is defined as: 


Q* = 


QL - Qh 


( 4 ) 


Qxx Qyy 

Q xx and Q yy are the reduced stiffnesses parallel and perpendicular to the fibers. Qgg is the 
transformed reduced stiffness corresponding to the angle 0 between the fiber orientations of 
the two plies adjacent to the delamination (Figure 2). The superscripts l and u refer to the 
lower and upper plies. For convenience, an expression for Q* in terms of the ply properties 
and the angle 0 is given in the Appendix. 


Dj e and Dq 9 are the flexural rigidities of the sections of the laminate above and below 
the delamination, respectively, in the direction of the fibers of the ply immediately below 
the delamination (Figure 2). These values can be found by considering the two sections as 
separate laminates and calculating their flexural rigidities according to [4] 


Dgg = 



( 5 ) 


where Qgg is the same as previously given, and z is measured from the midplane of the 
individual section, not of the entire laminate. 

The remaining three parameters er max , t m , and tf depend on the laminate as well as on 
the impact. cr max is the maximum out-of-plane tensile stress under the point of impact. t m is 
the time at which this stress occurs. Finally, tf is the duration of the stress. Wu and Springer 
[1] evaluated cr max , tf, and t m using a three-dimensional, transient, finite element method. 
Here, our analysis differs from that of Wu and Springer. We propose to approximate these 
three parameters in a simple manner, as described below. 
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The maximum stress cr ma x is taken to be proportional to the pressure P which arises at 
the flat interface of two elastic, semi-infinite cylinders when the two cylinders collide 


^max — k\P 

where k\ is a constant. The impact pressure P is given by [5] 

Pi V C,' 


(6) 


P = 


1 4 . BB. 

^ PcCc 


( 7 ) 


where V is the relative velocity between the two impacting surfaces, which, in our case, is 
the impact velocity, p is the density, and c is the dilatational wave speed in the material. 
The subscripts i and c refer to the impactor and the impacted cylinder, respectively. In our 
problem, the latter corresponds to the top ply of the composite plate. Thus, the wave speeds 
in the impactor and in the impacted composite may be calculated by the expressions 


c i — yfW/ Pi C c — \J Eyy/ ~Pc 


( 8 ) 


where Ei is the Young’s modulus of the impactor, and E yy is the transverse ply modulus of 
the composite. 

The time t m required to reach cr max is taken to be proportional to the time it takes for 
the first wave (caused by the impact) to travel in the composite from the impact surface to 
the back surface of the laminate. Thus, we write 


tm — t» 


( 9 ) 


and 


t w = h/c c (10) 

where ki is a proportionality constant and h is the thickness of the laminate. For plates 
made of sublaminates of different materials, this expression can be generalized to 


_ h& 
tw ~ 


(ii) 


where is the thickness of the sublaminate, is the wave speed in the material of 
that sublaminate, and n is the total number of sublaminates. 
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The duration of the stress if is approximated by the time during which the impactor 
is in contact with the impacted surface. 


tf tct 


( 12 ) 


This contact time is taken to be proportional to the contact time for the one- dimensional 
problem of a rigid body of mass m impacting a massless, elastic cylinder of length h and 
diameter 2 r (twice the nose radius of the actual impactor). 


tct = &3 t\D 

For the one-dimensional problem the contact time is [6] 


t\D 


ith / m \ 
c c \wr 2 hp c ) 


1/2 


(13) 


(14) 


For plates with sublaminates of different materials, the corresponding equation is 

1 1/2 


t\D = 


hW 


Tiif i r-^* 

r 2 l ( c 0>) 2 p ( j ) 


C J 


(15) 


where h^\ c£\ are the thickness, wave speed, and density of the sublaminate, 
respectively. 

By substituting Eqs. (6)-(14) into Eqs. (1) and (2) we obtain 

C 2 


to 

A and B are constants 


«■>*©* (Iff 

m _ B (M£] Ct m c ° (^) Cs C, “ 


V K Ic ) 


DJo j 


/ J, \ C*9 

B = k ° 7 (!f) C 2 
\«2/ 


(16a) 

(166) 

(17) 

(18) 


The constants A, B, C 2 -C 5 , and Cj-Cio must be obtained by fitting Eqs. (16a) and 
(16b) to experimentally determined delamination lengths and widths. Therefore, we per- 
formed a least squares fit of these equations to the data given in Reference 2 and summarized 
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in Tables 1 and 2. This procedure yielded the following values 

A = 180 B = 20 

C 2 = 3/4 C 7 = 4/3 

C 3 = 1/2 C 8 = 4/3 (19) 

C 4 = 1/4 C 9 = 1/3 

C 5 = 2/5 Cio = 0 

For these values, the expressions for the delamination length and width become 



As previously discussed, the multipliers (180 and 20) as well as the exponents are constants 
independent of material property, layup, and impact velocity. 

IV. CALCULATION PROCEDURE 

By using Eqs. (20) and (21) the delamination length and width can be estimated between 
any two plies in the plate. To perform the calculations we must first select the location (or 
locations) at which we wish to determine the dimensions of the delaminations. At each 
location of interest the delamination length and width are calculated by Eqs. (20) and (21). 

The parameters Kj c and are material properties and their values must be determined 
by static tests (Eq. 3). The value of Q* depends on the material, on the layup, and on the 
location under consideration, and can be calculated by Eq. (4). The flexural rigidities Djg 
and Dg 0 depend on the material and on the layup and can be calculated by Eq. (5). Note 
that the location of the delamination is taken into account through Dj e and Djg. The impact 
pressure P , the duration of the stress tip, and the time t w required to reach the maximum 
stress depend on the material, the layup, and on the impact, and are to be calculated by 
Eqs. (7), (10), and (14). 
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The calculations can be performed readily by a personal computer or even on a pro- 
grammable calculator. To perform the calculations we have written a computer code (des- 
ignated as “SIMPACT”). The input parameters required by the code are given in Table 3. 


V. DISCUSSION 

To illustrate the use of Eqs. (20) and (21), delamination lengths and widths were calcu- 
lated for the conditions which existed in the tests reported in Reference 2 and summarized 
in Tables 1 and 2. Essentially, in these tests, Fiberite T300/934 graphite epoxy plates (16 
plies thick, 4 in by 4 in) were impacted by 1/2 in diameter aluminum spheres [2]. The plates 
were inspected by an ultrasonic scanning technique (C-scan) before and after impact. After 
impact the plates were also cut along their axes of symmetry, and were examined by an 
optical microscope. The measured shapes of the delaminations are shown as solid lines in 
Figures 3-6. Delaminations lengths and widths were also calculated by the model (Eqs. 20 
and 21) using the material properties listed in Table 4. The results of the calculations are 
indicated by dotted lines. The comparisons in these figures show that the calculated and 
measured delamination lengths and widths agree closely. 

A more quantitative assessment of the accuracy of the model can be made by comparing 
directly calculated delamination lengths with measured delamination lengths and calculated 
delamination widths with measured delamination widths. Such comparisons are made in 
Figures 7 and 8 and in Tables 1 and 2. Ideally, all the data should fall on the 45 degree solid 
lines. Deviations from these fines may be due to errors either in measurements or in the 
model. It is encouraging that the measured and estimated delamination dimensions agree 
within 25 percent for both the length and the width. 

The foregoing results indicate that the present model can be used to estimate the lengths 
and widths of delaminations under a wide range of conditions. It is hoped that future tests 
will verify the model also to be valid under as yet untested conditions. 
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APPENDIX 

EXPRESSION OF THE REDUCED STIFFNESS MATRIX Q 


The components of the reduced stiffness matrix Q are [7] 


Qxx — 

Qxy — 

Qyy — 


Exx 

1 - VxyVyx 
^xy Eyy 
1 ~ VxyVyx 

hi 

1 “ V xy Vy x 


Qss — G X y 

Qee = Qxx cos 4 9 + 2{Q xy + 2 Q sa ) sin 2 0 cos 2 9 + Q yy sin 4 9 


(A.1) 

(A.2) 

(A3) 

(A4) 

(A5) 


where E xx is Young’s modulus in fiber direction, E yy is Young’s modulus transverse to fiber 
direction, v xy and v yx are Poisson’s ratios, and G xy is the shear modulus. 

By substituting the appropriate expressions for Q xx , Q yy , and Qg$, we obtain 

Q l xx - \Q U xx + 2 ( Qxy + 2 Qss) sin2 6 cos2 0 + Qyy sin4 *] , * ^ 

Q QL^QZ, {AM) 

The superscripts l and u refer to the lower and upper plies, respectively. 
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Table 1 

SUMMARY OF DELAMINATION LENGTHS OBTAINED BY C-SCAN 













Table 2 

SUMMARY OF DELAMINATION WIDTHS OBTAINED BY C-SCAN 
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Table 3 

SUMMARY OF THE INPUT AND OUTPUT PARAMETERS 
FROM THE “SIMP ACT” CODE 


Input Parameters 

1) Impactor 

a) Mass, m 

b) Nose radius, r 

c) Density, pi 

d) Young’s modulus, Ei 

e) Velocity, v 

2) Composite Plate 

a) Ply layup 

b) Material properties (see Table 4) 

Output Parameters 

1) Final delamination length 

2) Final delamination width 


Table 4 

MATERIAL PROPERTIES OF FIBERITE T300/934 
GRAPHITE/EPOXY (from Reference 2). 


Ply thickness, h 

6.250 x 10 -3 

in 

Density, p 

5.548 x 10“ 2 

lbm/in 3 

Longitudinal Young’s modulus, E xx 

2.109 x 10 7 

psi 

Transverse Young’s modulus, E yy 

1.450 x 10 6 

psi 

Shear modulus in x-y direction, G xy 

8.251 x 10 5 

psi 

Poisson’s ratio in the x-y direction, 0 v xy 

0.3 


Transverse tensile strength, Sjt 

8.006 x 10 3 

psi 

Critical stress intensity factor of mode I, 6 Kj c 

1.183 x 10 3 

lbf/in 15 


“From Reference 4. 
J From Reference 8. 




MASS: m 
RADIUS: r 
VELOCITY: V 




BOUNDARIES: 
SIMPLY SUPPORTED 
OR CLAMPED 


Figure 1. Description of the problem 
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Figure 2. Illustration of the delamination 
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Figure 5. Comparison between the measured and calculated delaminations. Solid lines: 

measured by C-Scan[2]. Dotted line: calculated delaminations (eqs 20 and 21). 
T300/934 graphite/epoxy plates ([O 4 /- 454 / 454 / 9 Q 4 ]) with clamped edges 
impacted by 1/2 inch diameter aluminum spheres. 
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Figure 6. Comparison between the measured and calculated delaminations. Solid lines: 

measured by C-Scan [2]. Dotted line: calculated delaminations (eqs 20 and 21). 
T300/934 graphite/epoxy plates ([O4/454/6O4/9Q4]) with clamped edges impacted 
by 1/2 inch diameter al umin um spheres. 
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Figure 7. Delamination lengths in 3 inch by inch T300/934 graphite/epoxy plates impacted 
by 1/2 inch diameter aluminum spheres. Dotted fine is the result of the least 
square fit of eq (20). SS: edges are simply supported; CC: edges are clamped. 




Figure 8. Delamination widths in 3 inch by inch T300/934 graphite/epoxy plates impacted 
by 1/2 inch diameter aluminum spheres. Dotted linpis the result of the least 
square fit of eq (21). SS: edges are simply supported; CC: edges are clamped. 



